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Abstract 

We study Euclidean Wilson loops at strong coupling using the AdS/CFT correspon¬ 
dence, where the problem is mapped to finding the area of minimal surfaces in Hyper¬ 
bolic space. We use a formalism introduced recently by Kruczenski to perturbatively 
compute the area corresponding to boundary contours which are deformations of the 
circle. Our perturbative expansion is carried to high orders compared with the wavy 
approximation and yields new analytic results. The regularized area is invariant 
under a one parameter family of continuous deformations of the boundary contour 
which are not related to the global symmetry of the problem. We show that this 
symmetry of the Wilson loops breaks at weak coupling at an a priori unexpected 
order in the perturbative expansion. We also study the corresponding Lax operator 
and algebraic curve for these solutions. 
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1 Introduction 

According to the AdS/CFT correspondence string theory on AdSs x S ' 5 background is dual to AT = 4 
super Yang-Mills theory in four space-time dimensions [I]. This theory is considered to be integrable 
in the planar limit, and over the years integrability based techniques where used and developed for 
making progress in solving the spectrum of the theory, see |2:j for review. Most of the study of 
integrability in this theory is related to the spectral problem. Much progress was also made in the 
context of the scattering amplitudes/lightlike Wilson loops duality. On the other hand, much less is 
known about integrability properties of Wilson loops, and in particular Euclidean Wilson loops. As 
is well known, at strong coupling the problem of computing the expectation value of a Wilson loop 
is translated to finding the area of a minimal surface in AdS space, ending on the AdS boundary, 
where the contour is defined by the Wilson loop [3|. In this paper we focus on the problem of finding 
this area in Euclidean AdS 3 or equivalently the three dimensional hyperbolic space, H 3 . 
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Over the years several approaches for studying the problem were proposed. The corresponding 
bosonic sigma model can be classically reduced to the cosh-Gordon model using Pohlmeyer reduction 
[41. The solution to the cosh-Gordon equation is known in terms of theta-functions [5j|7]. Thus, for 
any solution there corresponds a Riemann surface, or equivalently an (odd genus) algebraic curve 
which satisfies some special properties in order for the solution to be real. This algebraic curve was 
recently shown to be the same curve which is extracted from the Lax operator [8j. Unfortunately, 
given a Wilson loop contour, it is not known what is the corresponding algebraic curve. Similarly, 
given an algebraic curve, the properties of the Wilson loop, e.g. the area, contour etc. are a priori 
not known, moreover it is not even clear that the solution corresponds to a closed Wilson loop (for 
example, the solution may correspond to some minimal surface ending on an infinite self intersecting 
line, or ending on more than one contour, which might be unstable). These issues make it hard to 
use this formal solution for practical purposes as well as studying other properties of these objects. 

Another interesting approach was recently introduced by Kruczenski |9|. In this approach the 
worldsheet boundary is fixed to be the unit circle, and the problem is mapped to finding the correct 
parametrization of the contour in the conformal gauge. Once the correct parametrization is found, 
one has all the required ingredients to compute the area in terms of a contour integral around the 
unit circle on the worldsheet. The main drawback is that it is not known how to find the correct 
parametrization either analytically or numerically. 

A very different approach to the same problem in Minkowskian AdS 3 space was recently proposed 
in jlO|. This approach uses the known solution to the problem where the Wilson loop contour is 
given by a set of null lines found in |11| , where the nontrivial limit of infinite cusped segments is 
taken, so that the contour becomes smooth. The resulting area is given in terms of an integral TBA 
equation, which one should generally solve numerically. 

As can be understood form the discussion above, the solution to the problem is not expected 
to be simple by any means. More precisely, the only "simple" examples of known solutions are the 
circular Wilson loop and the infinite straight Wilson line. Beyond that, the simplest solutions are 
the "genus-one "[^solutions which correspond for example to the correlation function of two circular 
Wilson loops, the gg-potential and the infinite cusp, which are all given in terms of Jacobi elliptic 
functions. Other solutions should correspond to higher genus Riemanu theta-functions. In order to 
understand better the properties of these object it is desirable to have other simple examples at hand 
which can be studied analytically. One useful approach is the so called "wavy expansion 


which is an expansion in a deformation parameter around the circular or infinite line contours, 
keeping terms up to second order. However, this is not always satisfactory and some interesting 
features can be missed since to this order all the properties are fixed by conformal symmetry |12] , 

In this paper we use Kruczenski’s formalism introduced in |9| to study the problem perturba- 
tively around the circular contour. We also adapt this approach for studying perturbations around 
the infinite straight line contour. We show that quite easily we can go far beyond the wavy ap¬ 
proximation, and apply this approach to several simple perturbative deformations. Among these 
Wilson loop contours are some families of symmetric contours which include the ellipse and some 
less symmetric contours which include the limaqon. Besides for providing us analytic expansion for 
the area to very high orders in the expansion parameter, having this analytic data for these Wilson 
loops allows us to study some nontrivial properties of these objects. 

One aspect we shall study is the expansion of the Lax operator and the resulting algebraic curve, 
which hopefully can shed some light on the relation to the Riemann theta-functions. Another aspect 
is related to a curious property which was pointed out in j6j. There it was shown that there exists 
a one parameter family of deformations of the target space contour which leaves the area invariant. 


1 Here we refer to the genus of the corresponding Riemann surface. 
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The deformation depends on the spectral parameter A which must be a phase in this case, and is 
not related to any obvious global symmetry of the problem (i.e. global conformal symmetry). We 
shall call such deformations A -deformations throughout the paper. Because of the universality of 
the wavy correction due to conformal invariance 12 , it is expected that this deformation leaves 


the Wilson loop expectation value invariant to third order in the expansion parameter for any 
value of the coupling. However, this universality breaks at the fourth order in the perturbation 
of the circular Wilson loop (l4 and it is not obvious whether this invariance at strong coupling 
should survive quantum corrections, and in particular what happens at weak coupling. Using our 
approach we check how this deformation affects the one-loop weak coupling expectation value. We 
find a dependence on the deformation parameter in all of the examples we consider. Interestingly in 
all of the examples except one the dependence shows up starting at the 8 th order in the expansion 
parameter (in the limagon example it starts to show only at the 16th order), and the effect is very 
mild. Moreover, given an arbitrary contour, it is not trivial at all to find its A-deformed contour. 

Since there are no other analytic results available for the area of minimal surfaces ending on 
simple contours as we study in this paper (e.g. the ellipse), we can only compare our results to 
numerical data, as we do in the case of the ellipse where such numerical data in available |l5], and 
find agreement. Furthermore, our analytic results should be useful for testing new approaches to 
come for solving the general problem. 

Finally, although our main interest is the study of Wilson loops and their properties, let us 
mention that the same mathematical problem is of great interest also in the study of entanglement 
entropy of conformal field theories in 2 + 1 dimensions 1161 , where our approach and results can be 
also applied. 

The paper is organized as follows, in section [2] we start by giving a brief introduction to the 
formalism presented by Kruczenski in [9]. We provide only the relevant results which we shall 
use throughout the paper, and omit the derivation which can be found in |9|. In section [3] we 
perform a general wavy perturbation around the circular Wilson loop solution, to low orders in the 
wavy expansion, which should set the ground for the later analysis of concrete examples where the 
same procedure is applied to much higher orders. We then study some specific examples in detail 
in section [f] where the ellipse is set as our prime example for which we provide a more explicit 
analysis. Later, in section [5] we slightly modify the formalism such that we map the upper half 
plane to the minimal surface instead of mapping the unit disk. This is of course equivalent in 
principle, and we can study the same deformations in this gauge, however, it is more convenient for 
some deformations of the infinite straight line solution as we shall demonstrate. Finally, in section 
[ 6 ] we summarize and discuss our results. In appendix [A| we provide details relevant for the general 
solution of the area to fourth order in the expansion parameter. 


2 General setup 


In this section we briefly repeat the analysis presented by Kruczenski in |9|, where the reader is 
referred to for more details. We start with the bosonic string sigma model on H 3 , and the first step 
is to perform a Pohlmeyer reduction. We use the embedding coordinates X ij= q ^ in M 3,1 subject 
to X 2 = Xq — Xf — X '2 — X| = 1, where the SO(3, 1) ~ SL( 2,C) symmetry acts linearly. The 
Poincare coordinates are given by 


X -MY 


AT + iX 2 
X 0 - X 3 ’ 


X 0 - X 3 ' 


( 2 . 1 ) 


The worldsheet metric is taken to be Euclidean and the coordinates are given by z = a + ir, 
z = a — ir, and the corresponding derivatives are given by d = | (Jj^ — i-§y) and d = (^ + ij^). 
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The action in conformal gauge is given by 


S = 


1 

2 


dacLr {dX,,dX 11 - A(A>Y^ - 1)) , 


( 2 . 2 ) 


and it should be supplemented by the Virasoro constraints 

dX^JX 11 = 0 = dX^dX' 1 . (2.3) 

The equations of motion are given by 

ddX^ + = 0, (2.4) 

where consistency with the constraint X 2 = 1 yields A = BX^dX^. Following J6j we rewrite the 
embedding coordinates in terms of a matrix X = where a^ = {1 , cd} and a 1 are the Pauli 

matrices. In this way the constraints and equations of motion are given by 

X f = X, det (X) = 1, det (dX) = 0, det (<9X) = 0, BBX - J tr (BXa^dXa^) X = 0, (2.5) 

where we used the relation ^ tr (Acr^Bcr^). The first constraint is solved by introducing 

a new matrix A by X = AA^. The second constraint implies that det A = e 1 ^, but the phase 
can be cancelled by gauge fixing as we shall see immediately so we can assume det A = 1, so 
A G SL(2,<C). The equations are invariant under the global symmetry transformation X —> UXU t 
and A —» U A with U G SL( 2,C). Furthermore, there is a gauge symmetry A —>• A H(z,z) which 
leaves X unchanged if H(z,z) G 17(2). 

Next we define the traceless currents J = A _1 <9A and J = A~ l BA and rewrite the constraints 
and equations of motion in terms of the new variables. Using the equations of motion and constraints 
one can construct a flat connection given by 


A A) 


-\da fe a \ 

Xe a T^da J ’ 


7(A) 


Ida \e a \ 
—fe~ a da J ’ 


( 2 . 6 ) 


where a(z, z) and f(z ) satisfy the generalized cosh-Gordon equation 


dda(z,z) = e 2 “("’ s ) + \f(z)\ 2 e~ 2a ^\ 


(2.7) 


and where A G C is the spectral parameter. a(z, z) is a real function and f(z) is holomorphic 
where the solution is defined. For A = 1 we get back the currents defined by the string solution. 
Interestingly, in case where A is a phase this still corresponds to a string solution, generally defined 
by a different contour, with the same area for the minimal surface ending on the original contour. 
This deformation of the original contour is not related to any obvious global symmetry of the 
problem and we shall call such deformations A -deformations throughout the paper. Notice also 
that defining j (A) = J(\)dz + J{\)dz we have the relation 

j(A) t = -j(-A ) = (2.8) 

Next, following |9| we consider a minimal surface solution parameterized such that its contour 
in target space is mapped from the unit circle on the worldsheet. Using polar coordinates on the 
worldsheet z = re 1s , let us denote the contour by in target space by X(0) = X(0) + iY{9) where 
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X(9) and Y(9) are the Poincare coordinates at Z = 0, or equivalently at r = 1. Defining r = \Jl — £, 
close to the boundary the solution for a(z, z) is given by 

o) - - ine+£(i+ e)fa(o) +e 4 / 3 m + (2.9) 

where fa{9) are real functions, fixed by fa(9) and f(z). 

Very interestingly, assuming one knows the correct parametrization function F{9 ) for the bound¬ 
ary contour X(6 ), using the linear problem 

d'tj) = i/^J, d'tjj = iJ)J, (2.10) 


it is possible to derive the following relation 


1 


2 ^ 


{X\F(9)), 9} = - - I2fa{9) - 2\f(9)e zw + -f{9)e 


-2 id 


( 2 . 11 ) 


where {,} stands for the Schwarzian derivative. Y A=O (0) = X(9) is contour we start with, and 
for A being a phase, X x (9) is a deformed contour giving rise to a minimal surface with the same 
regularized area. Throughout the paper we shall parameterize A = e lv> with ip £ [0, 27t] . Equation 
(2.11) will be of prime importance in the current paper. The expression implies that if we know 


the correct parametrization in the gauge we are using, then in principal we have all the required 
information in order to reconstruct f(z ) and a(z,z). Moreover, using further manipulation (which 
can be found in |9|), Kruczenski presents the following expression for the regularized area 


M reg = —27T — 4 


drda\f\ 2 e 2a =—2 tt— 


- <b d9 


Re{Y, 9} - {w, 9} 


do In w 


( 2 . 12 ) 


where w(z) = j z \ffdz. and the integral is performed around the unit disk. Although this construc¬ 
tion is very appealing, at present it is not known how to find the correct parametrization either 
analytically or numerically. In the following we shall approach the problem perturbatively around 
the simplest available solutions, namely the circle and the infinite line, and show that this formalism 
enables us to perform such a perturbative expansion very efficiently without finding the minimal 
surface explicitly, and to reach very high orders in the perturbation expansion compared to the 
wavy approximation. 


3 Perturbations around the circular Wilson loop solution 


Using the formalism introduced in |9], we study perturbations of minimal surfaces around the 
circular Wilson loop solution. We start with a general analysis which yields a simple expression 
for the area to third order in the wavy expansion. This result was already obtained in 
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using 

different techniques, however we believe it is useful to re-derive it here in order to set the ground 
for higher order perturbations which we will compute later. We compare the result with the weak 
coupling computation and observe the expected universality up to this order. General higher order 
analysis can be carried in principle, however it is far more complicated since conformal invariance 
is less restrictive at higher orders. We give the general solution to the next order in appendix [A] 
however we do not find a simple general expression for the area in terms of a contour integral. In 
the next section we shall turn to some explicit examples and demonstrate the power of applying 
the formalism to obtain results to very high order in the perturbative expansion. 

We start by collecting the required properties of the circular contour solution which is the 
starting point of the later analysis. 
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3.1 The circular Wilson loop properties 

The circular Wilson loop is the simplest example which can be solved analytically, where the minimal 
surface is the semi-sphere. The (non-unique) A matrix which we will use is given by 


A c 


1 

\Jl — zz 



(3.1) 


From this we can extract a c = In f c (z ) = 0 and = 0, where the generalized cosh-Gordon 
equation is reduced to the simpler Liouville equation. The boundary contour in the conformal gauge 
is given by X(9) = e l9 and we have {X(9),9} = We can also solve explicitly for A C (A) 


a m _ 1 ( Q(A)^ + c(A) x(a(A) + c(X)z) \ 1 f 1 ±z\ 

y/1- zz\ a(X)z + c(A) y(«(A) + c(A)z) ) y/1 - zz \ A z 1 ) 


(3.2) 


where we picked a(l) = 0, c(l) = 1 and a(l) = l,c(l) = 0 for simplicity and concreteness. Thus, 
obviously in case where A is a phase, this corresponds to a rotation of the surface]^] Finally, the Lax 
operator is given b}0 


L 0 (z,z) 


A 

1 — zz 


1 + zz 2f A 
— 2Az —l — zz ) ' 


(3.3) 


Notice that the Lax operator is not unique, we choose it such that the resulting algebraic curve is 
given by 


y 2 = — det Lq = A 2 , (3.4) 

so that it coincides with the limit of the genus one solution up to a constant [7], namely 

y 2 — lim = A(A — a)(A-|—) 

a—>-oo d 

This is equivalent to the trivial curve y 1 = 1 found in 


3.2 General perturbation of the circular Wilson loop minimal surface 

In this section we find a simple general expression for the area of the wavy circle to third order, by 
adding a small perturbation to the circular contour. We compute a(z,z ) and F(9 ) to first order 
where F{9) is the parametrization function which will be properly defined below, and f(z) to second 
order, which fixes the area to third order. As we shall see, the reason that the third order area 
formula is fixed by the lower order functions is that f(z) vanishes to leading order for the circle. 
The fourth order result which is not fixed by conformal symmetry is much more complicated. In 
appendix |A| we given the general solution to a(z, z) and f(z) to the second and third orders which 
can be used to compute the area integral to fourth order, however we do not find a simple expression 
for this integral in terms of a contour integral. 

2 This should be expected since the deformation leaves the regularized area invariant, and there is only one contour 
with the regularized area equal to — 2-7T, namely the circle 

J Here we refer to the Lax operator of the Pohlmeyer reduced model which is related to the sigma model’s Lax 
operator by a simple gauge transformation |8|. 
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(3.5) 
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3.2.1 The general setup 


Our starting point is a given contour on the H 3 boundary which depends continuously on a parameter 
e, and reduces to the circle, X{9 ) = e 10 , when e = 0. Let us denote the contour by X(0) = 
e iF(0)+J2 n=1 £ n Gn(F(0)) w j iere the function F(9) is the correct parametrization function which is 
unknown. Notice that the given functions G n {6 ) are arbitrary complex periodic functions, G n (s ) = 
Sfcez G^e iks . Obviously, by these definitions F(9 ) should also depend on e and should reduce to 
F{9) = 9 when e = 0. Thus, we expand the parametrization function in powers of e as follows, 

m = 9 + Zn=ie n F n (e). 

Next, we wish to compute the e expansion of the minimal surface area ending on X(9), which 
to leading order is given by — 2ir, the area of the senri-sphere. In principle we need to find a(z, z) 


and f(z) which enter the area integral (2.12). We expand these functions as follows, 


1 ] 

a(z, z) = In --r = In -- 5 - + a n (z, z)e r ' 

1 — zz 1 — r z 1 

OO 


n =2 


(3.6) 


n =1 


Notice that there is no term linear in e in the a(z, z) expansion since it vanishes due to the boundary 
conditions. We decompose the perturbation to the real and imaginary parts as G n {9) = G r n {9 ) + 
iG l n (9), where G r n (9), G l n (9) G M. Our starting point is the Schwarzian derivative of the contour, 
satisfying 

Re{X(9),0} = ^-12/3 2 (0), 

Im {X(0),9} = -41m (e 2ie f(e i6 )). (3.7) 

The second relation implies that 

{X(9),9}), (3.8) 

where V projects on positive frequencies, that is V = ^(1 +T~L) with T~L being the Hilbert transform 
on the unit circle, normalized such that 

Ji [G(t)} (s) = p.v. J dtG(t ) cot > ( 3 - 9 ) 


for which we have T~L [e mi ] (s) = sign(n)e* 
cosh-Gordon equation (2.7). 


Finally, a(z, z) and f(z) are related by the generalized 


3.2.2 Finding f(z ) to second order 

Next, we are going to use these relations to find f(z) to second order, which is enough to fix the 
area to third order in e. To leading order we have 

Re{*(0), 9} = ^ + e£ 3 (*i(0) + G\(0)) + 0(e 2 ) = ^ 12^{0), 

Im{X(0),0} = -eCsG^e) + Q(e 2 ) = -41m (e 2W f(e ie )), 


(3.10) 







where we defined the operator £3 = + <%. Since /3 2 (9) = 0(e 2 ), we have 

Fl (0) = -G[(9) + F h (9), 


(3.11) 


where F h {9 ) = a + 6sin# + ccos# is the homogenous solution, which we shall dropQ so F\{9) 
—G\(9). From the second equation in (3.10) we read fi(z) using (3.8), so 


fi(z) = -e-™nCzG\me^-n. z , 

Using F\{9) given above, the e 2 term in Im{X(#),#} follows, 

Im {X(9),9} = -e£ 3 G[(9) - e 2 £ 3 ( G r 2 (9 ) - G^dgG^)) + 0(e 3 ), 
so f 2 (z) is given by 

f 2 (z) = l -e~ 2W V (£3 ( G r 2 (9 ) - G\{9)deG\{9))) 

This is all we need in order to fix the area to third order. 


(3.12) 

(3.13) 

(3.14) 


3.2.3 The area 

The regularized area expansion is given by 

Aeg = -27T - 2e 2 J d 2 zf 1 (z)f l (z)e ~ 2a ° ( *- 2) - 2e 3 J d 2 z (fi(z)f 2 (z) + h{z)f 2 {z)) e~ 2a ^ 

~ 2e 4 J d 2 z (fi(z)f 3 (z) + fi(z)f 3 (z) + f 2 (z)f 2 (z) - 2f 1 (z)f 1 (z)a 2 (z,z)) e ~ 2a ° (z ’ z) 

+ 0(e 5 ). (3.15) 


Let us define 

00 

Im{X(#), 9} = J2 e n AA(0), (3.16) 

n =0 


so that 

fn(z) = l -e- 2ie V(£ 3 S n (9))\ e ^-anz = -\ [zd 3 + 3d 2 ) (VS n (9)\ e ^-n nz ) = ~ (zd 3 + 3c> 2 ) G n {z), 

(3.17) 

where 


Gi(z) =iPG^(#)|^_ ilnz , 

G 2 (z) =V (G 2 (9) - G\(9)deG r 1 (9)) 

G 3 (z) =7>( (c r 3 - G\G 2 - G^Gi - \G^d{G\) 2 + l((G\) 2 Gl) 


+ £ 


-1 


1 


2iF £ 3 <a + kc 3 g\ - - (A) - 2g;g; £ 3 <a 


-1((c;') 3 + 3g;o;'c;") 


0 —>•—iln.2 


(3.18) 


4 One might want to keep the homogenous solution in some cases so that F( 0 ) =0 and F( 2tv) = 2-k. This can 

always be fixed by choosing the constant piece of F h appropriately, however £3 does not "see" this term so it doesn’t 
change the current analysis. 
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(See appendix | a| for definition of G 3 (z)). Now, as was shown to first order in |9], with these definition 
we generally have 




^d (fj{z) ( 2 zGi(z) + (1 - zz)(2dGi(z) + z( 1 - zz)d 2 Gi(z 


(3.19) 


so we can use Stokes theorem to convert the surface integrals over the disk into line integrals over 
the disk’s boundary, that is 

I d 2 zfi(z)f j (z)e~ 2ao ^ z ’^ =-j> dOe~ 2ie f j {e~ i0 )Gi{e ie ) = - % -j dOC^Gi, (3.20) 

or 

I d 2 z(f i (z)f j (z) + f j (z)fi(z))e- 2a °^ z ’^ =~lj> d0n{Si)C 3 Sj = -- j dOUiS^Si. (3.21) 
Before writing the area, let us mention that these integrals can be expressed also as 
J d 2 z(fi(z)fj(z ) + f j {z)f i {z))e~ 2ao ^ z] = j) ds <j> dt cot Si{t)£ 3 ,sSj(s ) 

= — ^ ds j) dt esc 4 ^ ^ Si(t)Sj(s). (3.22) 


where we wrote explicitly the Hilbert transform and integrated by parts to act on the kernel, instead 
of on the £) function^ Using the expressions above, we can write the regularized area as 


Ae g = -2n + e z -j) d9U{G\)C 3 G\ + e A i £ dOV.{G r l )£ 3 {G 2 - G\d 0 G\) + 0(e 4 ). 


(3.23) 


We did not hnd such a nice expression for the fourth order contribution, see appendix [Aj 

Finally, let us comment on the universality of this result. The weak coupling one-loop integral 
is proportional to 


tt- is / j/W ■*(») - l*(»)ll*(f)l 

W 1 = jdsjit (xW _ j(t))2 ■ 


(3.24) 


Throughout the paper we ignore the prefactoij®] which should appear in front of W\. It is 
straightforward to expand Wi to third order in e, and after some algebraic manipulations to arrive 
at 

Wi = -2vr 2 + 2vre 2 ^ j) dOU^)C 3 G\ + 2ire 3 i £ deH(G\)C 3 (G r 2 - G\d 0 G\) + 0(e 4 ). (3.25) 

Thus, the second and third order terms have the same form at weak and strong coupling so these 
terms are universal. The reason for that is conformal symmetry |l2[[l4|. It was shown explicitly 
in 1141 that the form of the fourth order terms is different. 


5 This form is more similar to the expressions one gets when computing the Wilson loop at one-loop at weak 
coupling. 

6 here A corresponds to the’t Hooft coupling and not the spectral parameter. 
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4 Applications 


In this section we study some specific examples in detail. We will examine several deformations of the 
circle with different symmetry properties. We will consider a family of contours which interpolate 
between the circle and hypocycloids. These hypocycloids have cusps, and we do not expect the 
expansion to be valid for values where the cusps become significant (and logarithmic divergences 
are expected), instead, the expansion should be valid when the contour is relatively smooth. One 
special case of these contours is the ellipse which we will consider in more details first. These 
contours have both reflection and discrete rotational symmetries, according to the number of cusps. 
We will also consider a different family of contours with the same symmetries as the previous family. 
After that we will consider a contour with less symmetries, the limagon which has one reflectional 
symmetry. At last, we will consider a contour with no discrete symmetries. 

Since the analysis for different contours is conceptually very much the same, we chose to present 
a more detailed derivation only for the ellipse. In the other case we only quote the final results. In 
all cases we present results to very high order in the expansion parameter relative to the ordinary 
wavy approximation. 


4.1 Ellipse 

We express the ellipse curve by 

X{8) = e iFW + iesmF(6), (4.1) 


where the parametrization function F{9) is a real monotonic function satisfying F{ 0) = 0 and 
F(2n) = 27r (for definiteness, although there is some freedom which allows more general boundary 
conditions), and e is a real small number, which will serve as an expansion parameter, e = 0 
corresponds to the circle where all the details of the solution for the minimal surface are known and 
provided in subsection 3.1 Next we are going to solve the problem order by order in e. First we 


take the parameterizations function to have the form 


F(9) = 0 + yF n (0)F 


(4.2) 


71=1 


with F n (0) satisfying F n ( 0) = F n (2ir) = 0. We shall rely heavily on the Schwarzain derivative 

erivative is given by 

((l + .)co,(F( g )) +i ,ln(F(»)))0 F '«>) a - + 2 n»)F"'(») 


relations (3.7), where in this case the Schwarzian derivative is given by 

3e(2+e) 


{X(9),d} = 


1 - 


2 F'{8f 


(4.3) 


where the primes represent a derivative with respect to 8. In the following we will need the real 
and imaginary parts of the Schwarzian derivative which are given by 


Re{X(0),0} 


(-> . 3e(2+e)(—(1+e ) 2 cos(F( 6 >)) 2 +sin(F( 6 >)) 2 ) \ 4 

^ ((l+e) 2 cos(F( 0 )) 2 +sin(.F( 6 O ) 2 ) 2 J W 

2 F'{8) 2 


3F"{8) 2 + 2F'{8)F'"{8) 


(4.4) 


Im{X(0),0} 


6 e(l + e)(2 + e) sin(2 F{8))F'(8) 2 
(2 + e(2 + e) + e(2 + e) cos(2 F{6))) 2 


(4.5) 


11 









From these equations we can read f(z) and /?2(0), 


h (e) = 


i i 


12 V 2 


-Re{X(0),0} , 


f(e id ) = - ~e~ 2ie V (lm{X(9), 0}). (4.6) 

Next we expand the real function a(r, 9) and the holomorphic function f{z) as, 

oo oo 

a{r,9) = ^Ta n (r,6)e n , f(z) = ^ f n (z)e n , (4.7) 

71=0 71=1 

where an = In 1 1 - = In . 

u 1—zz l—r z 

As for the general analysis, the procedure goes as follows. First we use the fact that given a 
solution to the problem to order e n , Im{ A, 9} is known to order e n+1 . This allows us to (easily) 
extract f(z ) to order e n+1 . Thus, we start with the circular Wilson loop solution and immediately 
read f\ = — Next we plug f(z) in the generalized cosh-Gordon equation and solve for a(r,9). 
We use the two implicit boundary conditions, a n >o( r = 1) = 0 and a n (r = 0) = finite. Then we 
expand a(r,9) around r = 1, extract 02(0) and plug in (4.6) to solve for F n (9). 

Due to the symmetries of the ellipse the expansions take the form 


1 = 0 + 5] F n {9)e n , 

71=1 

, F n (9) 

= Y^Fn,k sin 2 k9, 

k= 1 

(4.8) 

OO 

/0) = ^2fn(z)£ n -. 

71=1 

i fn(z ) 

71—1 

= '^2fn,kZ 2k , 
k= 0 

(4.9) 

OO 

= ^a n {z,z)e n , 

71=0 

a n (z,z) 

71—2 

= 5] a n,k( r ) COS 2k9, 
k= 0 

(4.10) 


a 


and it turns out that the expansion procedure is quite simple and can be carried to high orders. 
Moreover, we observe a pattern in the functional form of a(r, 9), 

i(6(n—fc) —ll+3(—l) n - fc ) 

a n ,k(r) = r 2k (l-r 2 ) 2 ^ a n ,k,ir 21 , (4.11) 

1=0 

which holds at least to the 18th order in e, and we suspect it might hold to any order. This reduces 
the problem to an algebraic one of finding constant coefficients. Finally, we find the following 
expression for the regularized area 


Tlreg(e) — 2-7T ^ 


37re 2 37re J 

+ 


2377re 4 | 1177re 5 64881vre 6 | 64443vre 7 14373577Tre 8 

4 320 + 160 89600 + 89600 20070400 

3584953vre 9 1103146882197re 10 220647325797re n 663090748836438l7re 12 

+ . .„. + 


5017600 154542080000 30908416000 

1106373532973931vre 13 40943000996733445243vre 14 


9281797324800000 


1546966220800000 57175871520768000000 

1952095942839819321vre 15 1 577506909298315380292446977re 16 


+ 

+ 

+ 27471862748298608640000000 

+ 0(e 19 ). 


2722660548608000000 219774901986388869120000000 

19736906966190071806502297Tre 17 801650044535506237372382994066703vre 18 


1115068403809909423032238080000000 


(4.12) 
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Figure 1: The area of the minimal surface ending on an ellipse contour parameterized by (cos(s), (1+e) sin(s)). 
The different lines correspond to different orders in the e expansion (from order e 2 up to order e 18 ), where 
it is seen how the lines converge. For e > 0 we plot the map of the solution for — 1 < e < 0 which behaves 
better at this region as explained. On top of the graph we added the points generated numerically in |15| . In 
the simulation in 


15 


, the authors parameterized the ellipse as 


R\ 


v 


= 1, and the blue points correspond 


to taking R\ = 1 and so our i-axis is e = (R 2 /R 1 ) — 1- The red dots correspond to taking R 1 = 2. As one 
can see the accuracy in these two cases in slightly different. Moreover, the numerics are limited by the IR 
cutoff (which was 0.03 in this case), so we do not expect a perfect agreement, although it is easy to see how 
our analytic prediction agrees with their data. We see that our approximation begins to break for e > 4, 
which correspond to R 2 /R 1 = 5 where the ellipse is already very different from the circle we started with. 


Notice the alternating sign in the expression which implies that A reg (e) gives a better approximation 
for — 1 < e < 0 rather than e > 0. Fortunately, due to conformal symmetry it is enough to consider 


we 


— 1 < e < 0 only, where the e > 0 values are related by e' = where — 1 < e < 0. In figure [l] 
plot the area as a power of e using different e power approximations from order e 2 to e 18 . For e > 0 
we plot the image value using (4.12) for e < 0. We also plot numerical data points on top of the 


functions using the data from 15 ^ It is also interesting to see how the parametrization function 
looks like, keeping in mind that it is enough to find this function in order to solve the problem. We 
plot the function in figure ([ 2 ]). 

As discussed in the previous sections, there exists a one parameter family of deformations of the 
ellipse contour, A-deformations, for which the area does is left invariant. Given the analysis above 
it is relatively easy to find these contours X * to the same order in e, by solving 


{X\6), 6} = Re{A A=1 (0), 6) - 2A f(0)e zw + -f{0)e 


-2 i6 


(4.13) 


where f(6) = f(r = 1,6), with the "initial" condition X A=O (0) = X{9). Since A must be a phase, 
we define A = with p G [0, 27t]. We plot the deformed contour for ip = 0, 7t/47t/2, 37t/4, n in 
figure ([3]). Notice that for the circle this deformation amounts to rotations and does not give new 
solutions. 

'I am grateful to P. Fonda and E. Tonni for sharing this numerical data. 
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Figure 2: (a) The parameterization as a function of s for e = —1/2 (which is equivalent to e = 1) up to order 
e 17 (b) F{6) — 0 as a function of 8 for e = —1/2 (which is equivalent to e = 1) for all orders up to order e 1 '. 







Figure 3: A deformations for the e = —1/2 (or equivalently e = 1) ellipse. We use A = e lLp and in (a) p = 0, 
(b) tp = j (c) p = | (d) ip = ^ and in (e) <p = tt. The computation is carried to order 0(e 13 ). 


4.1.1 Weak coupling and A-deformation 

One can speculate whether the A-deformation of the Wilson loop contour leaves it invariant beyond 
the strong coupling limit, a' quantum corrections are hard to compute, however, the one loop weak 
coupling quantum correction is easier to evaluate. Next, we compute the one loop weak coupling 
correction to the for the A-deformed ellipse X x (0), by expanding in powers of e. We take A to be a 
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Figure 4: The relative deviation of the one-loop weak coupling expectation value of the A-deformed ellipse 
with <p = 7 t/2 compared to the undeformed ellipse, using equation (4.141. As can be seen, the difference is 
very small. Notice that the expectation value of the deformed Wilson loop is smaller than the one of the 
ellipse for any p. 


phase, A = e*^, and denote by a subscript the p> dependence. We find 


W ltV ,= 


ds 


x(t)-x(s) - |a(g)||a(t)| 
(x(s) — x(t)) 2 


2 3vr 2 e 2 37T 2 e 3 10l7r 2 e 4 53vr 2 e 5 17597r 2 e 6 1805vr 2 e 7 

I I I 

2 ~2 64 ~32 1024 1024 

TT 2 (8227199 + 7776cos(2(^))e 8 ( tt 2 (2078399 + 7776 cos(2 ip))e 9 
4587520 + 1146880 

vr 2 (4632279678637 + 40879629488 cos(2v9))e 10 
2543321088000 

vr 2 (927152406637+ 15013543088 cos(2p))e n 
+ 508664217600 

tt 2 (740574919262795417 + 18893380267080208 cos(2<^>) + 30047328000000 cos(4(^>)e 12 

407338305454080000 

tt 2 (368445596660191451 + 13339409812485424 cos(2<^>) + 90141984000000 cos(4v?))e 13 
+ 203669152727040000 

+ 0(e 14 ), (4.14) 


The invariance to third order is expected since the coefficients to that order are universal due 
to conformal symmetry. However starting from order e 4 the coefficients should not be fixed by 
conformal symmetry alone. Quite surprisingly, the first appearance of A dependence shows at the 
8 th order. Thus, eventually the symmetry is broken by quantum corrections. Interestingly, the 
relative difference is extremely small as can be seen is figure [4] where we show the relative difference 
between W\^=q and W\^ =v j 2 - where it is maximal. 


4.1.2 Lax operator 

Finally, it is interesting to study the Lax operator and the resulting algebraic curve. It is known that 
minimal surface solutions in H 3 are related to theta functions defined on a hyperelliptic Riemann 
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surface of odd genus [5]-[8]. However, the relation between a given Riemann surface and a contour 
on the AdS boundary is poorly understood. Given a generic algebraic curve, it is not known 
a priori what is the corresponding extremal surface (or even if a closed smooth nonintersecting 
extremal surface exists) and vise verse. For example, it is not known what is the Riemann surface 
corresponding to the minimal surface ending on the ellipse, what is the genus and even if it is finite. 

Next we construct a Lax operator for the ellipse by solving the equation 


dL+[J,L\ = 0, (4.15) 

order by order in e, starting with the circular solution, where the Lax operator is given by 


Lq{z,z) = 


A 

1 — zz 


1 + zz 
—2A z —l — zz 


2 " 


(4.16) 


The Lax operator is not uniqiifQ but the resulting algebraic curves for different Lax operators should 
be related by a bi-rational transformation. At each order in e we find some integration constants 
which we can set as we wish. We choose these constants such that we arrive at the simplest result. 

The explicit form of the Lax operator for higher orders in e is not very illuminating, however, 
simple algebraic curve is found by taking its determinant. Each time when we increase the power 
of e the resulting Lax operator develops an extra A ” 1 factor in the L \2 component (similar to the 
one in Lq if we ignore the overall A pre-factor), which increases the order of the pole by 1 (a similar 
observation was made in [l 8 ]| where the general wavy approximation was considered.). However, 
the Lax operator remains a rational function of the spectral parameter which implies it is a good 
Moreover, we can set the integration constants and an overall constant factor such 


operator 

that 
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y 2 = — det L = 1 + 0(e 16 ). (4.17) 

Thus, we do not observe any interesting structure which can point to a relation with the Riemann 
surface corresponding to the exact solution. We comment more on this issue in the discussion 
section. We found the same result, namely a trivial curve, in all the examples we have looked into, 
thus we shall not repeat this fact in the following subsections. 


4.2 Hypocycloids 


In this section we apply the above procedure to a family of symmetric contours which interpolates 
between the circle and hypocycloids, which are n-cusped symmetric contours (see figure [5j, given 
by 


X{0) = e iF{e) + i ee ^-n/ 2 )F{e) sin (p F ( 0 )) } 


(4.18) 


where we recover the ellipse for n = 2. Accordingly, the contours have a symmetry under n — 1 
discrete rotations and n reflections. The cusps appear when e = —2 /n, and due to conformal 
symmetry if is enough to consider — 2 /n < e < 0 (for non-intersecting contours). 

Next, we give the regularized area of the minimal surfaces and the one loop weak coupling 
expectation value for the A-deformed contours for n = 3,4, 5 (the n = 2 case is the ellipse which 
we already considered). We plot the results for the area, correct parametrization function and 
A-deformations in figures §0 and [ 6 ] respectively. 

8 For example, any power of the Lax operator yields another Lax operator as well as any of their linear combinations. 
In the case of the circle Lq = 1 which is a "trivial" Lax operator, more generally we could have started with aLo + b. 
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(a) 






(d) 


Figure 5: Hypocycloids to circle deformation. Here we plot some examples of (4.18) for n = 3,4, 5,6, 
for different values of — - < e < 0. For clarity, the figures are scaled to avoid intersections. 


A n=3 

^reg 


A n=A = 
■'•‘■reg 


= — 2n — 3ire z + Sire 3 — 


12937T6 4 873vre 5 1386017vre 6 1948851vre 7 


+ 


+ 


280 ' 140 156800 ' 156800 

88672365813vre 8 31600967429vre 9 36254366537313759vre 10 in . 

+ _ _ __ + 0(e 10 ), 


5022617600 

157T6 2 157T6 3 

— 27T---1- 


1255654400 


1005528043520000 


(4.19) 


2109vre 4 1689vre 5 

+ 


13473123vre 6 24353319vre 7 

+ 


2 2 112 56 224224 224224 

63592443793533vre 8 299021266130137re 9 29627266774400450319vre 10 in . 

+---+ 0(e 10 ), 


305249584640 


76312396160 


39477684581823488 


4 71=5 
^reg 


— 27T — 157Tt + 157re 3 — 

+ 


176 ' 88 

7582306709205vre 8 4338088174725vre 9 


10335vre 4 9015vre 5 43653157re 6 9245505vre 7 

+ -XX-T^XXX- + 


15488 


15488 


W?= 3 = - 2 7T Z - 6? tV + 6vrV - , ( 

vr 2 (2133269689 + 3648094 cos(2 </?))e 8 
50226176 

7r 2 (764606393 + 3648094 cos(2<^))e 9 


5075603456 


2 2 


+ 0(e w ). 


-2,3 


+ 


1268900864 

8l7r 2 e 4 57vr 2 e 5 83vr 2 e 6 2 377r 2 e 7 


(4.20) 


(4.21) 


+ 


8 


+ 


12556544 


+ 0(e lu ), 


(4.22) 


Tirr) _ 4 2 2 9 2 '( 13337T 2 e 4 10937T 2 e 5 71179vr 2 e 6 129937vr 2 e 7 

W\ l - A = - 2vr 2 - 15vrV + 15vrV-—-1-—-—-h 


32 


16 


512 


512 


7r 2 (3837945912753851 + 7831242589272 cos(2v?))e 8 


+ 

rn =5 _ o—2 


7814389366784 

7T 2 (1816117657651195 + 7831242589272 cos(2 (p))e 9 


+ 0(e 10 ), 


1953597341696 

-2,3 ion—2,4 


w?=° = - 2 7T Z - 30vr 2 e 2 + 307r 2 e 3 - 1307r 2 e 4 + 230vr 2 e 5 - 

5vr 2 (15805224661 + 34793379 cos^t^e 8 9 
22658944 h ^ 


2585vr 2 e 6 5515vr 2 e 7 

-:-h- : - 


(4.23) 


(4.24) 


As is easily seen, in all of these examples the week coupling expansion of the A-deformed contours 
is dependent of the deformation, where the dependence starts at the eighth order. 
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Figure 6: A-deformations. We plot the A-deformations for ip = 0, f, for the interpolating contours 

between the circle and hypocycloids (ordered from left to right). We used e = —0.35, e = —0.25 and e = —0.2 
for the n = 3,4, 5 cases respectively. 


We can also find the area for general p, however in this case we could compute it only up to the 
fifth order in e, 


A reg = — 27r- n (n 2 — l) 7re 2 H —n (n 2 — l) 7re 3 

8 8 

n (56 + 8n - 307 n 2 + 8n 3 + 338n 4 - 160n 5 + 57n 6 ) vre 4 
512 (4n 2 - 1) 

n (24 + 8n - 147?r 2 + 8n 3 + 210n 4 - 160n 5 + 57n 6 ) vre 5 
256(4n 2 - 1) 


+ 0(e 6 ). 


(4.25) 


The reason is that for higher orders the number of terms that a(z, z) contains is proportional to n, 
which complicates the analysis. However, for a fixed n this is not a problem. 


4.3 Symmetric wavy contours 

Another family of symmetric contours we study is given by 

X(0) = e iF(0)+*dn(pFm t (4.26) 


see figure [9j 

as in the previous example, we can study it for different specific values of p. Here we consider 
two examples, p = 2 where quite easily we can get to very high orders in the expansion, and see how 
the weak coupling integral depends on the A-deformations in very high orders. The second example 
is p = 13, which is a bit harder, but one can still see the same dependence on A to the first orders 
in the expansion. Next we give the results 
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(a) 


(b) 


(c) 


Figure 7: Parametrization. We plot the correct parametrization for the interpolating contours between 
the circle and hypocycloids. (a) n = 3, e = —0.35, (b) n = 4, e = —0.25, (c) n = 5 e = —0.2. 





Figure 8: Regularized area. We plot the regularized area for the interpolating contours between the circle 
and hypocycloids. The different lines correspond to different orders of the e expansion, from e 2 to e 10 for 
n = 3,4 and e 9 for n = 5. We plot for — ^ < e < 0, where the cups form at e = — ^ and the area should 
diverge, (a) n = 3, (b) n = 4, (c) n = 5. 






(a) 


(b) 


(c) 


(d) 


Figure 9: Wavy contours to circle deformation. Here we plot some examples of (4.261 for p = 2, 3,4,5, 
for different values of 0 < e. For clarity, the figures are scaled to avoid intersections. 
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Figure 10: Regularized area. We plot the regularized area for the interpolating contours between the 
circle and the wavy contours. The different lines correspond to different orders of the e expansion, from e 2 
to e 18 for p = 2 and e 10 for p = 13. (a) p = 2, (b) p = 13, 





Figure 11: A-deformations. We plot the A-deformations for <p = 0, J,§,^p 
e = 0.5. 


7r for the p = 2 case with 




^reg 2 = — 27T - 37T6 2 + 


93vre 4 501437re 6 510139vre 8 65754318359vre 10 

+ 


+ 


+ 


20 4200 14400 582120000 

1195458440855851vre 12 61047851487256409vre 14 


3178375200000 47344547250000 

45707069078388982419341507vre 16 52566325973037148254959546391187vre 18 


10124976097716480000000 
+ 0(e 20 ), 


3273637646841985463040000000 


(4.27) 


= - 2vr - 1092vre 2 + 


16609327re 4 3887594024353vre 6 67968797564585251lTre 8 


25 570000 

2652706006393624451200787779Tre 10 12 . 

-b 0(e 12 ). 

24329522800000000 v ’ 


+ 


821712000 


(4.28) 
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Wf~ 2 = - 2vr 2 - 6vr 2 e 2 + 


31vr 2 e 4 9857r 2 e 6 vr 2 (10132841 - 11016 cos(2( / j))e 8 


+ 


161280 


+ 


+ 


4 48 

vr 2 (-507075437357 + 519046632 cos(2</?))e 10 
2483712000 

vr 2 (-205243491659573369 + 193885849040584 003(2^) + 6619276121760 cos(4<^))e 


12 


vr 2 e 14 


1628955430502400000 
+ 207824411640450672 cos(4</?) ) + 


298343485440000 

- 3872331127046431211507 + 3825980252246834960 cos(2v?) 


7r 2 e 16 


584860157767581696000000 
- 4894520980993299391859564659 + 5286197530092023750882320 cos(2v?) 

+ 304417890855492833823264 cos(4^) + 4656911757360684307200cos(6v?)^ + 0(e 18 ), 

27534945l7r 2 e 6 


Wf~ L6 = - 2tt 2 - 2184ttV + 1079267r 2 e 4 - ^ 

13vr 2 (64529639666077152031 + 16900003369957344 cos(2 tp))e' 


+ 


585469800000 

137T 2 (17188687999825875665056111457 + 5500782068869693320779168 cos(299))e 10 


1160518237560000000 


+ 0(e 12 ). 


(4.29) 


The p = 2 example is quite interesting since we could compute W\u> to higher orders than in 
the other examples above and see how the <p dependence changes as we increase the order of the 
expansion. A plot of the regularized area is given in figure [lO] In figure 11 we plot the contour for 
different values of the phase A = e llf> for p = 2, for p = 13 it is very hard to notice the change by 
the deformation for the values of e for which our approximation is valid, thus we do not present the 
corresponding figures. 

4.4 Limagon 

In this case the contour is given by (see figure [l2| (a)) 

X(6) = e iF{6) + ee 2iF{9) . (4.30) 

As can be seen, this contour has only one reflectional symmetry. The regularized area of the minimal 
surface is given by 

8977re 8 8347re 10 87254l7re 12 8200237re 14 

~20 

i nQQnon^onQQi n~,20 /moQ/i'70o r 7n„,22 

+ 0(e 24 ). 
(4.31) 


A re g = — 2-7T — 37re 4 — 127re 6 — 


5 1400 350 

69449016l7re 16 329884137re 18 193892 972033 1 97re 20 371402 3 472279vre 22 


78400 


980 


150920000 


7546000 


The one loop weak coupling expectation value for the A-deformed contour is given by 


W lv = - 2„ 2 - ft** - 247r 2 e 6 - - 3467r 2 e 10 - 21083lHe " - 20117 ^ 

,v 4 16 4 


2 ,14 


3tt 2 (115201973 + 2592 cos(2vj))e 16 tt 2 (83104799 + 7776 cos(2 ¥ 7))e 
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+ C>(e 20 ). (4.32) 
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Figure 12: The limagon. (a) The limagon for different values of e between e = 0 (the circle) and e = 1/2 
(the cardioid). For clarity, the different contours in the figure are scaled so that they do not overlap, (b) 
The minimal surface area calculated using (4.311 by including terms up to e" terms from n = 2 (upper blue 
line) to n = 22 (lower purple line), (c) The parametrization for e = 0.3125, corresponding to the 6 th inner 
contour in (a). 



In figure 12 (b) and (c) , we plot the regularized area and the correct parametrization function 
respectively. In figure 13 we plot the contour for different values of the phase A = e %v . The relative 
difference between and Wi tV>= o is very small. For example, for e = 0.3125 the maximal relative 

difference is ^’A'r Wl ’ , f=^/ 2 ~ 379 x iq -9 . 

w i,ip=o 

4.5 Asymmetric contour 

In this case we choose the simple asymmetric contour (see figure [l4| (a)) 

X(s) = e im + e ( e 2iF{e) + ie 3iF(0) ) . (4.33) 

The regularized area of the minimal surface is given by 


A-eg = — 27T — 37re z — 


19177re 4 350823vre 6 2475105369vre 8 


+ 0(e 10 ). 


20 200 78400 

The one loop weak coupling expectation value for the A-deformed contour is given by 

773vr 2 e 4 573597 r 2 e 6 tt 2 (1182155647 + 62208 cos(2(^))e 8 


IV[ w = — 2i — 67 T 2 e 2 — 


16 


17920 


(4.34) 


+ 0(e 10 ). 

(4.35) 
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Figure 14: The asymmetric contour, (a) The asymmetric contour for different values of e between e = 0 
(the circle) and e = 0.2. For clarity, the different contours in the figure are scaled so that they do not overlap, 
(b) The minimal surface area calculated using (4.341 by including terms up to e n terms from n = 2 (upper 
blue line) to n = 8 (lower green line), (c) The parametrization for e = 0.08, corresponding to the 3rd inner 
contour in (a). 







(a) (b) (c) (d) (e) 

Figure 15: A-deformations of the asymmetric contour for e = 0.07. (a) ip = 0, (b) <p = f, (c) tp = §, 

( d ) V = X’ ( e ) V = 7r - 


In figures 14 (b) and (c), we plot the regularized area and the correct parametrization function 
respectively. In figure 15 we plot the contour for different values of the phase A = e ltp . The relative 
difference between Wi tV> and W\^=q is very small. For example, for e = 0.07 the maximal relative 
difference is Wl ' v ^ v Wl ’ , f=^/ 2 ^ 197 x iq -9 . 


5 Perturbations around the infinite straight line 

In some cases it might be more convenient to map the upper half plane (UHP) to the target space 
Wilson loop. For instance, this could be the case for the infinite straight line solution, which is 
conformally equivalent to the circle. Furthermore, for the circular solution this is the gauge used 
to write the minimal surface solution in terms of theta-functions J 6 J. In this section we adapt the 
formulation introduced in |9] to map the solution from the UHP instead of the unit disk, and provide 
an explicit example of perturbation of the straight line into a sinusoidal curve. 

For the infinite straight line solution we take the worldsheet to be the UHP, where it is more 
convenient to use the a and r coordinates. The boundary curve is defined on the real axis, that is 
at t = 0. Using the fact that d a = d + 8 and d T = i(d — 5), the generalized cosh-Gordon equation 
is given by 

\{dil + 9 T 2 )a(a, r) = e 2a ^ + \f(a + ir)\ 2 e~ 2a ^\ (5.1) 
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Similarly, the flat-connection components are given by 

0 —a i 1 


MX) = 


2 o T a 

\e a - fe~ a 


f e 


MX) = 


— \d a °t 

i(Xe a + fe~ a ) 


iif MM 1 ' <5 ' 2) 


Notice that a(d, r) has to diverge at r = 0, so expanding it around that point using the EOM gives 

a ~ — ln(2r) + r 2 7 2 (cr) + r 4 74 ( 0 -) + 0(r 5 ) (5.3) 

where 74 (d) = ^ (l 6 |/(d )| 2 + 272 (d) — 72 (d)). In order to have a single boundary curve we should 
require e 2 Q ( <T > 7 Y 0 ) < 1001 (remember that this factor is related to the induced matric). 

Defining H(a,r) = and expanding d r H close to the boundary we find H(a, 0) = iX, 

which implies that the boundary curve is given by X = 

From the linear problem, by defining x = rir ' l lx ■ we get ~x” + ^( cr ) r )x = 0- Close to the 

boundary we have 


V(v) - ^72(d) + A/(d) - i/(d). 


(5.4) 


Next, we can relate the potential to the Schwarzian derivative by {X(a),a} = —2 V(<j,t = 0), so 


(X A (d), d} = - 372 (d) + 2 A/(d) - 2 ^/(d), 


(5.5) 


thus, 


Re{X A (d),d} = - 372 (d) 


1 


Im{X (d),d} = -2*(A/(d) - -/(d)) = —4Im(A/(d)). 


Now, 


(5.6) 

(5.7) 


/(d) = i(l + 77)Im(/(d)) = -- (1 + ‘R)Im{X(d),d}, 
where 7i [/(d)] (d 7 ) = — ip.v. f^° da represents the Hilbert transform along the line. 

5.1 Perturbation around the straight WL minimal surface 

Using the equations introduced above we can repeat the wavy analysis around the straight line and 
present high order results. We start by defining the straight line solution, and continue with an 
example of sinusoidal perturbation)^] 


5.1.1 Straight Wilson line 

The simplest solution is given for the infinite straight line where a(a, r) = — ln(2r), f(a + it) = 0. 
The Target space solution is given by 

/ — 2CT-fr (T — IT \ 

A = f \/ 2t / k A \/? r / K 


X = 


k(ct 2 +t 2 ) 

T 

Act 


At 

1 


iX 


( 5- 8) 

\ T KT / \ yj 2t K \/2 TK / 

so the target space Poincare coordinates are given by Z = kt and X = K,X~ 1 a (we kept the 
deformation parameter which in this case rotates the line in the X — Y plane). 


9 As a check, we also repeated the ellipse analysis using this parametrization for the circle, and got the same results 
we found in the previous section. The explicit solution looks quite different, but is still simple enough in order to 
carry the analysis to high orders. 
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5.1.2 Sinusoidal Wilson line 


Using the equations above we add a perturbation to the straight line. Generally, one can add a 
perturbabtion and solve the equations of motion and Virasoro constraints order by order. However, 
the procedure is complicated and it is hard to get to high orders. As in the previous section, the 
equations above simplify the analysis considerably, where one does not have to find the complete 
surface at each order. Here we show how it works for the simple case of the sinusoidal line where it 
turns out to be possible to compute analytically very high order corrections. The curve we would 
like to study is X(a) = a + iesincr. The procedure we use goes as follows. We stat with a general 
ansatz for a(a, r), f(a + ir) and the parametrization F(a) expanded in powers of e. At each order, 
we first "read" / from the imaginary part of the Schwarzian derivative, then we plug it in the 
generalized cosh-Gordon equation solving for a and use the boundary conditions at r = 0 and 
t —> oo to fix the integration constants. Finally we compare with the real part of the Schwarzian 
derivative and fix F. 

We find that for this case the expansion is given by 

oo n— 1 

a(a, r) = — In 2r 2 + ^ e 2n ^ cos{2ko)a 2n ^(r), 

n=l k=0 

oo n 

f(a + ir) = Y j e 2n ~ l £ f 2n . lt2k . u 

n—1 k=1 

oo n 

F(a) = cr + ^2 e2n sin ( 2ka ) F 2n,2k- (5-9) 

n=1 k =1 


a 2 n,k( T ) are simple functions involving T-exponentials times polynomials, that can be found ana¬ 
lytically. We can also write a(a, r) in a more explicit form as 


oo n— 1 

E' 2n 

n=1 k =0 


a(a,r) = - In 2r 2 + e 2n ^ cos(2ka) ^ e 2sT P 2s -k(r), 


(5.10) 


s=k 


where P n (T) is a polynomial of degree n. 

After finding these functions, the area can be easily integrated to give 


reg 
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662835985271c 10 750989224006549c 12 
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110592 18874368 10871635968000 15028949562163200 


14 


155198331721663659127c 

3665727786540072960000 


16 


+ 


132953910399497800665047c 

3619644351509283471360000 


+ 0(e 18 ). 


(5.11) 


where A is half the length of the worldsheet, o E [—A, A], This expression converge quite fast for 
e < 1. We can interpulate a function for the coefficients in order to estimate the area close to 


e = 1_. The results are presented in figure 16 We can also find the A deformed curve order by 


order, see figure 17 The deformed curve is very close to the original one, the difference shows first 
only in the e 3 correction term, while before that the deformation amounts to a translation of the 
curve in the X direction. 


6 Summary and discussion 

Minimal surfaces area expansion In this paper we studied minimal surface solutions in H3 
ending on various contours, by expanding perturbatively around the circular solution. We used the 
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Figure 16: The area of the minimal surface ending on an sinusoidal contour parameterized by 
(F(a), esinF(a)). The different lines correspond to different orders in the e expansion (from order e 2 up 
to order e 16 ), where it is seen how the lines converge. The dashed line is plotted using interpolation of the 
expansion coefficients up to e 1000 . As one can see, the wavy line approximation (blue line) is quite good for 
e < 0.4. 



Figure 17: The blue line represents the original sinusoidal curve for up to order e 17 for e = 0.9. The purple 
curve represents the resulting deformed curve with the same area, for A = e* 71 ’/ 2 (which gives the maximal 
deformation). In the plot we shifted the deformed line by a factor of ip = tr/2 in order to ease the comparison 
between the curves. 


formalism introduced recently in |9| to obtain high order results for the area compared to the wavy 
approximation, without finding the explicit minimal surface solution. We applied this approach 
to several target-space contours with different discrete symmetry properties. We explored several 
explicit examples including contours which interpolate between the circle and hypocycloids (which 
include the ellipse), some other wavy symmetric contours, the limaqon curve and a contour with no 
discrete symmetries. In the case of the ellipse which was treated in more details, we compared our 
results with the numerical data given in [l5] and found very good agreement in the range where 
our expansion is supposed to hold. In all cases we found quite easily expressions for the area, far 
beyond the wavy approximation (at least to order e 8 and up to order e 22 in one of the examples). 

Having such high order expansions, it is natural to wonder if one can express the area in a closed 
form or find a general expression to all orders for specific examples such as the ellipse. In the case 
of the ellipse we conjectured the functional dependence to all orders in the expansion which might 
help in this direction. 

A very interesting generalization of our analysis would be to study the problem in higher di¬ 
mensional space-times such as AdS 3 x S 1 or AdSd >3 where very little is known compared to the 

AdS 3 . 


26 








A-deformations As explained in ( 6 ], 9j, for any given boundary contour X(0) we can define a 
one-parameter family of contours X(6,ip) with ip E [0, 2n] where X(0, 0) = X(0,2ir) = X(9), 
such that the regularized area of the minimal surfaces corresponding to X(6,cp) is the same. The 
parameter (p is related to the spectral parameter by A = e zv> , and we referred to such deformations 
as A-deformations. The origin of this symmetry of the Wilson loops at strong coupling is not 
clear since it is not related to the global conformal symmetry of the problem, furthermore it is 
not clear if the symmetry survives quantum corrections. Moreover, given a contour X(6), it is not 
known what is the shape of X{Q,ip) without having the analytic solution at hand. We computed 
X (6, p) perturbatively in e and checked the dependence of the one-loop weak coupling Wilson loop 
expectation value on ip. Because of conformal symmetry one should not expect any dependence to 
order e 3 , however, at order e 4 there is no obvious reason for the weak coupling result to be invariant 
under the deformation. For all the examples we checked, it turned out that the ip dependence starts 
to show only at order e 8 except for the Limagon where it starts at order e 16 . Thus, we show that the 
symmetry does not survive quantum corrections in general. This fact also demonstrates the need 
for computing high order corrections. These results definitely could not be found using the wavy 
approximation (that is to order e 2 ) and not even by going to order e 7 . 

The ip dependence we found raises some natural questions. Although we checked only a number 
of simple cases, the fact that we find a dependence only at the eighth order (or higher) is quite 
suggestive, and it should be very interesting to understand if this is true for general deformations of 
the circle, and if this is the case what is the symmetry that protects from lower order corrections. 
Another natural question is how does the Wilson loop depends on ip at higher loop orders at weak 
coupling, and for the quantum corrections at strong coupling; when does the symmetry breaks and 
at which order in e? It is also quite interesting why the dependence showed only at the 16th order 
for the limagon which is a multiple of 8 which appears in all the other examples. 


Algebraic curve and theta-functions As discussed in ( 6 j, minimal surface solutions in H 3 
correspond to Riemann surfaces of odd genus defined by an algebraic curve y 2 ( A) = 0, which are 
invariant under A —> — 1 /A (in order for the solution to be real). Non-trivial smooth solutions which 
end on one closed compact contour (different from the circular Wilson loop) start to appear from 
genus three. Given a string solution, the algebraic curve can be found by constructing the Lax 
operator, which is a rational function of the spectral parameter, and taking its determinant |l9, 8 ]. 
We applied this procedure for our examples by computing the Lax operator perturbatively in e. 

In principle, one could imagine it is possible to deform existing cuts and add/remove cuts in 
order to continuously deform the target space contour until the desired contour is reached, and 
particularly to end up with the circular contour. However, this is not the picture our construction 
yields (see also 18 for analysis at leading order). We started with the trivial curve y 2 ( A) = 1 
corresponding to the circular contour, and at each order we were able to construct a rational Lax 
operator such that y 2 { A) = 1 remains. Although it was checked only up to a finite order in e for 
each example, it would be quite surprising if this result will suddenly change at higher orders. It 
would be very interesting to understand why do the algebraic curves we hud are trivial and whether 
a nontrivial algebraic curve can be extracted from our analysis. 

This may be due to special properties of our starting point, the circular solution, compared 
to the higher genus solutions. First we note that for the circular solution f(z) = 0, while the 
theta-functions are found when the f(z) term in the generalized cosh-Gordon equation is set to 1, 
after a change of coordinates followed by a gauge transformation. However, when f(z) vanishes 
such a transformation cannot be carried and the analysis becomes more subtle. In contrast, by 
construction our expansions have a well defined limit to the circular solution. 

In terms of the Riemann surface, the circular solution is a limiting case of the genus one solution 
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where the length of the two cuts shrinks to zero, namely, it is a genus zero solution with two poles. 
The limit of adding a pair of new cuts on top of this solution is not obvious in terms of the theta 
function solution and it should be interesting to see if such a procedure can generate new minimal 
surface solutions. 

Moreover, the worldsheet picture of the contour (i.e. zeros of the theta functions) is very 
different compared to the known higher genus solutions. For the circle, the worldsheet contour is 
a finite piece of an infinite straight line which is periodic in the target space. On the other hand, 
in the case of higher genus solutions the worldsheet contour is a closed loop which is mapped to a 
closed contour in target space. Thus, it is not clear if and how the worldsheet straight line can be 
continuously deformed to such closed loops. Notice that this is different from the genus one case 
which corresponds to two concentric circles contour in target space. There, it is possible to make the 
target space contours wavy by adding two cuts [7] or more and the straight lines on the worldsheet 
which are mapped to the two concentric circles become wavy. Therefore, it would be interesting to 
check if our analysis can be repeated for these solutions, and if it yields an algebraic curve which 
is consistent with the known exact results. In this case there are several complications, first the 
equations involve the Jacobi elliptic functions which are much harder to solve. Moreover, now one 
has to take care of two independent contours and allow for f(z) to be meromorphic inside the unit 
disk since the solution is defined in the annulus. For example in the case of two concentric circles 
f(z ) oc z~ 2 , although if one maps the solution from the UHP (or more precisely from the finite 
strip) f(z) is a constant as expected since the solution coincides with the theta-functions solution 
(see e.g. |20|). 
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A Higher order general analysis 

In this appendix we give the general solution for a(z, z) to second order and f(z) to third order. 
These expressions are much more complicated than the lower terms in the expansion, since conformal 
symmetry is less restrictive. For this purpose we have to solve the generalized cosh-Gordon equation 
to order e 2 

4(1 - r 2 ) 2 |/i(z)[ 2 + 8 ^ 2 ^2)2 “ ~2 0~ 2 d 2 + rd r + <9j) a 2 (r,d) = 0. (A.l) 

We can solve this equation by expanding fi(z) = /1 nZ n , where the coefficients f\ n are given 

by 

fl,n = — + l)( ra + 2)( n + 3)G'i n+2 , (A.2) 
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where G\ n is a coefficient defined by G£(s) = n eins ■ Plugging the f\{z) expansion, we 

find the general solution 


«2 (r, 9) = (r 2 - 1) el<>{n q) h,nh,q 

71=0 <J=0 


r^+n+q (4(3 + q) + n{ 4 + g))r 2+n+ ' i ' 

(3 + n)(3 + q ) (2 + n)(3 + n)(2 + g)(3 + q) 


+ 


2 (n + q + 4) 


Min(n,g) 


(1 + n)(2 + n)(3 + n)(l + q)( 2 + g)(3 + g) 


(q + n + 1 — 2k)r 


n-\-q—2k 


(A.3) 


/c=0 


From here we extract /? 2 ($) which is needed in order to find the parametrization, and is given by 
expanding around r = 1, 


h{9) AE/.,/./ 1 -"' 

n=0q =0 


__ n > a 

3(l+n)(2+n)(3+n)’ n - q 

_ ■»~2( 1 +g) _ > 

3(l+g)(2+ g )(3+g)> 9 - 71 


Notice that when n = q the expressions coincide. Let us further define 


n=0 g=0 
oo oo 

71=0 q=0 


q ~ 2(1 + n) 

3(1 + n)(2 + n)(3 + n) 

n — 2(1 + q) 

3(1 + g)(2 + g)(3 + g) 


(A.4) 


(A.5) 


We can use these expressions and the Hilbert transform in order to express j3 2 (9) without using the 
conditional expression, in order to finally express it in terms of the perturbation, 


fo{9) = \(B + + B~) + l -U(B + - B~) 


(A.6) 


Notice that there is no double counting of the zero modes since the Hilbert transform kills them. 
Next we insert f\ )U in terms of the G \ n 's using the relation (A.2) which gives 


1 

B + (0) =I5EE G W K!, '^ e " <<,+! ”(« - 2 < 1 + »))(! + «)(2 + «)(3 + 9), 

n =0 9=0 
1 oo oo 

fi'(»)= s EE G'l.n+2 el( ' n * 2)9 Gl q+2 e- < q+ V e (n - 2(1 + g))(l + n)(2 + n)(3 + n). (A.7) 

7i=0 q =0 


Since B ( 9 ) = H + ($) we shall concentrate only on B + {9). These expressions are related to the 
derivatives of G as follows, 

B+(9) = J- (2 dGC 3 G + GC 3 8G ) , (A.8) 

where here G = G\ with the zeroth and first modes removed. The relation with the Schwarzian 
derivative implies 

W) = -^3 (ci + F 2 - ±d(G\) 2 + \d{G\) 2 ^ + ^ (2 dG\C 3 G\ + G[C 3 dG[) . (A.9) 
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(A. 10) 


Defining F 2 = K — (G\ — \d(G\) 2 + \d(G\) 2 ) (where K = K(9)) we get 

C 3 K = -12p 2 (e) + * (2 dG\C 3 G\ + G\C 3 dG \). 

After some manipulations we find 

K = ^£ 3 1 (H [2dG\C 3 H[ + G\C ? /)H \]), 


(A.ll) 


where H\(0) = U[G\(9')\(9), i.e. G\ = g[(9) + g[(9) so H{ = g[(9) - g[(6), with g\(9) = 
Y^=i9i,n eind ■ Notice that C^ 1 e m6 = —i n (i-n : 2 ) , so for our periodic functions 


C^f(9) = e 


ie 


d9'e~ 2id ' 


d9"e ie " I d9'"f{9'") 


= / dBf(0) ~ o ( e 


-10 


d9e ie f(9) + e lti I d9e~ lti f(9) 


je 


-w. 


Notice also that £3 and its inverse commute with ~H. Let us now define 

OO 

lm{X(9),9} = Y,^^S n (9), 

n =0 

and plug F\(9) and £ 2 ( 0 ) in Im{ A ($),$}, which to third order gives 

£.3^3 ( 0 ) = - 

+2 k'C 3 G\ + KC 3 G( - i {(Gi) 2 - 2 G\G r ”^ C 3 G\ - \ ((G^f + 3 ) 


(A.12) 


(A.13) 


£3 ( Gl - G\G r 2 - G\Gi - \G r ;d(G\f + \((G\) 2 G^) 


+ 0(e 4 ), 


(A.14) 


now f 3 follows from (3.8). 


A.l Alternative derivation for the ct 2 (z,z ) solution 

We notice there exist an explicit particular solution for a 2 (z, z ) in terms of the perturbation 
a ( 2 p) (z, z) = WW - where W ( Zj ^ = 1 f z ( z2 _ _ 2 qv> 

( -L ZZ ) Z \ 


(A.15) 


where G(z) is related to f(z) and is dehned in (3.17). The homogenous solution which should be 
added (the one which is regular at r = 0 , as is) is given by 

a^\z, z) = | + ~~ (T(z) + T(z )) + zdT(z) + zdT(z), (A.16) 


where T(z) is an arbitrary holomorphic function and T(z) is its complex conjugate. Expanding the 
particular solution around r = 1 yields 


(p) _, zG\dG\ + zG\dG\ 

2 ~ 2(1 - r) 


1 


r—>l 


2(1 ~r) 


^2(n + q)^ n -^ e G n d q . 


n,q 


On the other hand, 


» 


a, ~ 


T + r 


1 — r 


— T 

1 — r 


AkQrp | ^—ikOrp 

e lk~\- e Ik 


~ ' k 

which fixes T in terms of G using the boundary conditions. 


(A.17) 


(A.18) 
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